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CONDITIONING  ELECTRON  BEANS  IN  THE  ION-FOCUSED  REGIME 


I.  INTRODUCTION 


The  ion-focused  regime  (IFR)  offers  several  advantages  for  propagating 

1-3 

and  transporting  intense,  relativistic  electron  beams.  Chief  among 

these  is  the  ability  to  propagate  in  a  self-pinched  state  free  from 

interactions  with  destabilizing  plasma  electrons.  Such  propagation  occurs 

when  the  beam  travels  through  a  plasma  less  dense  than  the  beam.  The 

radial  electric  field  of  the  beam  expels  the  plasma  electrons,  leaving  the 

heavy  positive  ions  to  confine  the  beam  electrostatically.  The  plasma  can 

2-10 

be  created  using  the  beam  to  ionize  a  background  gas,  or  using  a 
laser  or  a  pre-existing  lov-energy  electron  beam  for  ionization. 

The  first  approach  is  termed  "passive"  IFR,  while  the  latter  two  approaches 
using  external  ionization  are  termed  "active"  IFR.  In  either  case,  the  gas 
must  be  sufficiently  rarefied  that  the  plasma  electrons  can  escape  without 

4 

colliding  with  and  further  ionizing  the  gas  molecules. 

Several  IFR  applications  have  already  been  employed.  One  application 
is  to  transport  intense  beams  over  long  distances  without  external  magnets. 
IFR  transport  over  distances  exceeding  80  m  on  laser-produced  channels  has 
been  demonstrated  on  the  Advanced  Test  Accelerator  (ATA)  at  Lawrence 

1  1  A 

Livermore  National  Laboratory  ’  and  on  the  EPOCH  facility  at  Sandia 
National  Laboratories.^^  Both  of  these  facilities  used  an  excimer  laser  to 
ionize  a  trace  organic  gas.  A  second  application  is  suppression  of  beam 
instabilities  in  an  induction  accelerator.  IFR  transport  on  ATA  allowed 
the  solenoidal  magnetic  field  to  be  turned  off  throughout  most  of  the 
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accelerator,  resulting  in  a  factor  of  3-5  increase  in  the  amount  of  beam 

power  that  could  be  delivered  without  disruption  from  the  beam-breakup 
13  14 

instability.  '  Other  IFR  applications  include  beam  extraction  from 

18  19  19  22  23 

foilless  diodes,  ’  beam  confinement  *  and  extraction  in  cyclic 

24 

accelerators,  and  final  focus  in  linear  colliders. 

In  this  paper,  we  consider  yet  another  IFR  application:  namely,  to 

prepare  or  "condition"  a  high-current  beam  prior  to  its  injection  into  a 

dense  gas.  Experiments  have  demonstrated  that  an  unconditioned  beam 

straight  from  an  accelerator  or  diode  is  disrupted  by  the  resistive  hose 

25  26 

instability  after  propagating  only  a  short  distance  in  the  gas.  ’  One 

strategy  for  combating  this  instability  is  to  center  the  beam  and  reduce 

the  transverse  displacements  that  seed  the  instability.  '  A  second 

strategy  is  to  tailor  the  beam  emittance  so  that  it  decreases  from  head  to 

tail.  Such  tailoring  detunes  the  instability  and  can  reduce  hose  growth 
8  23  3A 

substantially.  ’  ~  In  general,  both  centering  and  tailoring  are  needed 

to  propagate  intense  beams  over  long  distances  in  dense  gas. 

The  IFR  mode  can  be  used  for  both  centering  and  tailoring.  Preionized 

12-14 

IFR  channels  are  especially  effective  at  beam  centering.  Centering 

can  also  be  achieved  using  passive  IFR  in  a  narrow  pipe,  as  was  first 

demonstrated  on  the  Experimental  Test  Accelerator  (ETA)  at  Lawrence 

Livermore  National  Laboratory^  and  later  on  the  SuperlBEX  facility  at  the 

35 

Naval  Research  Laboratory.  Passive  IFR  in  a  large  pipe  does  little 

centering  but  tapers  the  beam  radius.  The  radius  tapjr  can  be  converted  to 

emittance  tailoring  by  passing  the  beam  through  a  final,  thick  scattering 

36 

foil.  This  approach  was  first  used  on  ETA  and  has  been  the  primary 

g 

emittance  tailoring  technique  for  propagation  experiments  on  ATA  and 
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SuperlBEX.  ’  Passive  IFR  tailoring  has  also  been  employed  at  Sandia 
National  Laboratories.^^ 

Although  IFR  conditioning  has  been  shovn  to  improve  beam  propagation 

dramatically,  detailed  IFR  analysis  has  been  lacking.  In  this  paper,  ve 

develop  and  refine  analytic  models  and  design  formulas  for  assessing  the 

utility  and  limitations  of  the  ion-focused  regime  in  conditioning  intense 

37 

laboratory  beams.  Support  for  this  work  is  given  in  another  paper  that 

presents  fully  electromagnetic,  particle  simulations  of  IFR  conditioning. 

Analytical  models  of  IFR  propagation  have  been  developed  by  a  number  of 
3  A  23  38  AO 

investigators.  »  *  »  “  Here  we  concentrate  on  beam  centering  and 

radius  tailoring.  IFR  radius  tailoring  is  modeled  by  extending  the 

3 

equilibrium  analysis  of  Briggs  to  include:  emittance  growth  in  the  IFR, 

the  role  of  the  foils  used  at  the  transitions  into  and  out  of  the  IFR,  and 

better  treatment  of  the  non-equilibrium  beam  head.  Beam  centering  is 

analyzed  in  terms  of  the  deflection  forces  produced  on  a  beam  by  an  ion 

channel  and  the  container  walls,  and  the  growth  in  emittance  caused  by 

these  forces.  Ve  conclude  that  under  proper  conditions,  active  IFR  cells 

are  effective  at  centering  while  passive  IFR  cells  are  effective  at 

tailoring.  Limitations  from  ion-channel  collapse,  the  ion  hose 

instability,  ion-impact  ionization,  and  magnetic  trapping  of  the  plasma 

electrons  are  discussed  as  well.  Not  discussed  are  effects  like  gas 

38 

scattering  and  Ohmic  erosion  of  the  beam  head  which  are  usually 
unimportant  for  laboratory  applications. 
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II.  RADIUS  TAPERING  IN  A  PASSIVE  IFR  CELL 


A.  EquilibriuB  Model 

3 

Briggs  developed  a  simple  model  for  predicting  the  equilibrium  radius 

a^(T)  of  a  beam  that  creates  its  own  IFR  channel  through  beam  impact 

ionization.  Here  ve  extend  that  model  by  allowing  the  beam  current, 

emittance,  and  energy  to  vary  with  time  t  into  the  beam  pulse.  Ve  make  the 

following  assumptions:  (i)  the  beam  is  highly  relativistic  with  a  Lorentz 

factor  Y  »  1;  (ii)  the  plasma  electrons  are  created  solely  through  beam 

impact  ionization  and  are  ejected  instantaneously;  (iii)  the  plasma  ions 

are  immobile;  (iv)  the  radial  profile  of  the  beam  density  is  flat  out  to 

radius  a^;  and  (v)  a^  decreases  with  t.  The  assumption  that  the  beam 

profile  is  flat-topped  (and  therefore  self-similar)  is  not  rigorously 

justified  but  greatly  simplifies  the  analysis  without  seriously 

compromising  its  validity.  For  consistency,  we  retain  this  assumption 

throughout  most  of  the  analysis.  The  assumption  that  the  plasma  electrons 

are  collisionless  has  been  studied  by  Yu^^  and  discussed  by  Buchanan.^®  In 

general,  the  assumption  is  valid  at  gas  densities  below  a  critical  value, 

38 

given  in  room- temperature  air  by  P^(torr)  =  (Ij^/lO  kA)(l  cm/aj^).  Effects 

from  ion  motion  and  magnetic  trapping  of  plasma  electrons  are  discussed  in 

Sec.  VI.  Finally,  because  we  are  mainly  interested  in  experimental 

applications  where  the  gas  density  and  propagation  distance  are  small,  we 

38 

neglect  gas  scattering  and  Ohmic  erosion. 

Given  our  model  assumptions,  plasma  ions  are  created  according  to 


an. 
_ 1 

ax 
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"b  ’ 


(1) 


where  n.  and  n,  are  the  beam  and  ion  densities,  respectively,  and  v.  is  the 
0  1  1 
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rate  at  which  the  beam  ionizes  the  gas.  The  rate  is  proportional  to  the 
density  and  atomic  number  of  the  gas,  but  depends  only  weakly  on  the  beam 
energy  for  y  »  1;  in  room- temperature  air,  =  P(torr)  ns~^.  For  a  flat- 

topped  collapsing  becun,  we  can  rewrite  Eq.  (1)  in  terms  of  the  volumetric 
charge-neutralization  fraction,  f.  =  n^/nj^,  as 


rr  •  I?  ■  '’i  -  hh 


■  "i  -  hh  i"(V“b')  • 

2 

Here  -  ecna^  n^^  is  the  beam  current.  The  plasma  electrons  are  neglected 
because  they  have  been  driven  outside  the  beam. 

Ve  eliminate  the  dependence  on  a^  by  combining  the  definition  of 
normalized  rms  beam  emittance. 


c  =  yau 
n  b 


(3) 


ymc 


with  the  radial  pinch  condition  for  the  average  beam  temperature, 


1,42 


Ti  =  efjljj/2c.  (4) 

Equation  (3)  assumes  a  flat-topped  beam  and  neglects  radial  variations  in 

_2 

Y,  while  Eq.  (4)  neglects  self-forces  to  order  y  <<  fj  and  balances  beam 

pressure  against  the  ion  pinch  force.  For  a  flat-topped  beam,  the  rms 

radius  a  and  half-current  radius  a.  are  related  by  a^„^  =  a, zo  = 
rms  1/2  rms  1// 

ajj/'fl.  Although  most  results  are  given  in  terms  of  a^,  some  results  are 

37  8  9 

given  for  convenience  in  terms  of  Simulation^  and  experimental  ' 

results  commonly  use 


5 


Combining  Eqs.  (3)  and  (4)  produces 


V-^^i^b 


(5) 
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where  ■  me  /e  =  17  kA  is  termed  the  Alfven  constant.  (The  Alfven 

/  O 

current  is  traditionally  defined^  as  times  the  Lorentz  factor  PyO 
Inserting  Eq.  (5)  into  Eq.  (2)  yields 


This  can  be  immediately  integrated  to  yield 


fi(T) 


[■ 

2  'Ty  I 


dx' 


4y  Ii 


b  0 


(6) 


(7) 


Inserting  Eq.  (7)  into  force-balance  equation  (5)  produces 


V.  -Ty 

4~  c~ 
n 


>(7  I 


b 

A 


(8) 


Observe  that  no  assumption  has  been  made  regarding  how  or  when  the 

emittance  e  was  generated, 
n 

Result  (8)  is  based  on  the  assumption  of  equilibrium,  which  is  usually 
justified  only  if  the  IFR  cell  is  long, 


z  >  0.5  Xj  .  JrV^i'b  ’  <’) 

where  X_  is  a  betatron  oscillation  wavelength.  Early  in  the  beam  head, 

P 

condition  (9)  is  not  satisfied  because  fjljj  "^0.  As  a  consequence,  result 
(8)  overestimates  the  degree  of  radius  tailoring  actually  achievable.  A 
correction  is  given  in  the  next  section. 
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The  equations  for  and  a^^  simplify  if  and  y  are  constant.  In  that 
case,  Eq.  (7)  reduces  to 


f.(T)  =  0.5  v.Qjj(t)/Ijj(t)  , 


(10a) 


while  Eq.  (8)  reduces  to 


2e 


%<■')  • 


JrVjQ^<T)/I, 


Here  the  time-dependent  beam  charge  is  given  by 


Qjj(r)  =  JdT'  Ijj(T') 

0 


(10b) 


(11) 


Result  (10b)  demonstrates  that  radius  tailoring,  <  0,  occurs  as 

increases  monotonically  with  time  x.  For  example,  if  rises  as  x"’,  the 

charge  and  radius  vary  as  Qj^  «  x"''’^  and  aj^  “  ^-(m+l)/2^  respectively. 

3  -1/2 

Briggs'  result,  a^^  «  x  ,  is  recovered  by  setting  m  =  0. 

Another  useful  example  is  a  beam  current  that  rises  and  falls  as  a 
half-sine  wave. 


Ijj(x)  =  sin(iix/Xp) 


(12a) 


where  is  the  peak  current  and  x^  is  the  pulse  duration.  Then, 


ajj('c)  ®  [  l-cos(nx/Xp)  I 


-1/2 


(12b) 


The  beam  thus  pinches  rapidly  at  early  x  but  levels  off  at  late  x.  In  this 
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case,  the  neutralization  fraction  becomes  large,  »  1,  as  t  -»  Once 

f.  exceed  unity,  the  plasma  electrons  are  no  longer  ejected  but  are  pulled 
back  by  the  plasma  ions.  Nevertheless,  the  analysis  remains  approximately 
valid  because  the  plasma  electrons  ejected  earlier  oscillate  at  large 
radius  about  the  ion  column,  rather  than  re-concentrating  within  it.  As  a 
consequence,  the  pinch  force  on  the  beam  is  still  largely  controlled  by  the 
ion  line  density  filb/c. 


B.  Beam  Head  Dynamics 

The  preceding  model  is  valid  for  those  portions  of  the  beam  that  have 

attained  equilibrium,  but  it  is  not  valid  for  the  weakly-pinched  beam  head. 

A  simple  correction  is  to  assume  free  expansion  in  the  head  and  equilibrium 

in  the  body.  This  composite  model  provides  a  good  estimate  for  the  maximum 

head-to-tail  radius  variation  produced  from  a  passive  IFR  cell.  A  more 

44 

accurate  approach  would  be  to  solve  the  Lee-Cooper  envelope  equation. 

Ve  compute  the  radius  of  the  freely  expanding  beam  head  by  assuming 

that  the  beam  is  injected  with  a  Maxwellian  thermal  distribution  determined 

by  the  initial  emittance  and  rms  radius  rj^(O).  The  unpinched  head  then 

44 

evolves  to  a  Gaussian  profile  expanding  at  a  constant  rate: 


e  z 
o 


rrh(O)  ' 


(13) 


Here  z  is  the  propagation  distance,  and  we  have  assumed  rj^(z)  >>  rj^(O).  Ue 
multiply  the  rms  radius  r^^  by  'T?  to  compute  the  corresponding  edge  radius 
for  the  flat-topped  profile  assumed  in  the  model.  Free-expansion  and  wall 
scrape-off  thus  restrict  the  edge  radius  to 
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(14) 


a  Min{b,  2ze^/Ya^}  , 

where  b  is  the  pipe  radius  and  a^  «=  'f2rj^(0)  is  the  initial  edge  radius. 

Restriction  (lA)  shows  that  the  radius  of  the  early  beam  head  varies 
with  z  but  not  T,  assuming  constant  E^/ya^.  In  thi„  region,  the  solution 
to  Eq.  (2)  is  fj(T)  -  twice  the  value  given  in  Eq.  (10a)  for 

a  collapsing  beam.  This  additional  growth  raises  f.  in  the  equilibrium 
region,  so  that  Eq.  (10a)  should  be  replaced  with 


-  fi(^i)V^i> 


=  0.5  v.(Qjj(t)+Qjj(x.)] 


(15a) 


in  the  equilibrium  zone,  T  >  To  incorporate  this  change,  we  insert  Eq. 

(15a)  into  Eq.  (5)  to  obtain  a  first-order  revision  to  Eq.  (10b): 


2e 


aj,(x) 


YVi[Qb(x) 


Qb(^i)]/lA 


2e 


(x)/I^l  ^  irB^/zy 


(15b) 


Here  we  have  ignored  beam  scrape-off  and  have  normalized  ajj(O)  to  the  tree- 
expansion  limit,  2ze^/Ya^.  This  approximation  applies  to  both  the  pinched 
and  unpinched  regions  of  the  beam.  The  transition  timv:;  x^  to  equilibrium 
is  given  by 


ViQb(Ti)  = 


(ao/z) 


(16) 
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C.  Eaittance  Growth  fron  Anhanonic  Pinching 

Result  (8)  demonstrates  that  temporal  variations  in  alter  the 

radius  tailoring  obtained  from  a  passive  IFR  cell.  In  this  section  ve  show 

that  the  IFR  cell  itself  changes  the  emittance.  Ve  conclude  from  this 

calculation  that  a  passive  IFR  cell  produces  radius  tailoring  only  if  the 

beam  is  injected  warm,  Tj^  >  0.  If  the  beam  enters  cold,  it  exits  with 

little  radius  tailoring  but  with  the  emittance  inverse  tailored,  3c^/3t  > 

0.  As  discussed  in  the  appendix,  such  a  beam  is  predicted  to  be  highly 

hose  unstable  when  injected  into  dense  gas. 

A  simple  estimate  of  emittance  growth  for  an  initially  cold  beam  is  as 
45 

follows.  If  the  IFR  pinch  force  is  harmonic,  each  beam  electron 
undergoes  sinusoidal  motion  so  that  its  time-averaged  rms  radius  contracts 
by  1/^.  At  the  same  time,  the  rms  radial  velocity  of  each  electron  rises 
from  zero  to  r^cOp/'TJ,  where  r^  is  the  injection  radius  and  tOg  is  the 
oscillation  frequency.  For  a  uniform  ion  channel, 


“8  ■ 


(17) 


In  general,  the  beam  and  ion  channel  are  not  flat-topped,  so  the  pinch 
force  is  not  harronic.  The  electrons  then  do  not  oscillate  at  the  same 
frequency,  and  they  become  out  of  phase  with  one  another.  As  orbit 
coherence  is  lost,  the  beam  electrons  thermalize  to  a  steady  rms  radius  and 
transverse  thermal  velocity  given,  according  to  the  above  analysis,  by 


and  r^^co^Z-JT,  respectively,  where  r^^^  is  the  rms  radius  at  inject! 


bo 


ion 


and  is  an  average  oscillation  frequency.  This  suggests  that  the 
P 

normalized  emittance  of  a  beam  in  a  fixed,  anharmonic  centering  force 
evolves  via  phase-mixing  to  a  value  no  smaller  than 


^(min) 

n 


^boV'^- 


(18) 
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This  result  may  also  be  derived  by  relating  the  final  transverse  beam 
temperature  to  the  work  done  by  the  pinch  force  in  contracting  the  beam.^^ 
The  latter  approach  is  particularly  useful  if  the  pinch  force  is  highly 
enharmonic. 

Result  (18)  and  the  discussion  preceding  it  indicate  that  a  passive  IFR 
cell  does  not  produce  radius  tailoring  in  an  initially  cold  beam.  Instead, 
the  beam  contracts  uniformly  by  l/-(2,  and  its  emittance  rises  with  t  as 
fjljj  rises.  To  produce  radius  tailoring,  the  beam  should  be  preheated  so 
that  the  initial  emittance  e  satisfies 


e 

o 


> 


e(min) 


(19a) 


To  produce  head  expansion,  the  emittance  should  additionally  satisfy 


c  >  ya  . 

o  o 


(19b) 


According  to  Eq.  (15b),  radius  tailoring  will  then  occur  up  to  a  time  > 
Tj,  where  is  defined  by 

ViQb(%)  =  I8(e^/Ya^)^  -  (a^/z)^]  .  (20) 

For  X  >  x^,  the  radius  remains  roughly  co;.stant  at  a^^  =  a^/JT,  and  the 
emittance  grows  as  the  ion  pinch  strength  increases  according  to 

^i^b  =  ^i<^e>^b<^e>  ^  -  Qb(%>5 
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III.  ENTRANCE  AND  EXIT  FOILS 


A.  Foil  Heating 

IFR  cells  usually  require  entrance  and  exit  foils  to  maintain  the 

necessary  gas  pressure.  These  foils  scatter  the  beam  and  increase  its 

47 

angular  divergence  by 

<eS  =  <e.S  +  .  (22) 

2 

Here  <0^  >  characterizes  the  initial  angular  spread  of  the  beam,  t^  is  the 
foil  thickness,  and  is  a  scattering  coefficient  that  depends  on  the  foil 
material.  Foil  scattering  thus  raises  the  normalized  mean-squared  beam 
emittance,  *  0.5  by 

+  0.5  ,  (23) 

where  €.  is  the  initial  emittance.  Because  k^  is  nearly  independent  of  the 
beam  energy  or  current,  foil  scattering  can  serve  as  a  convenient  means  for 
producing  constant  emittance  in  a  pencil  beam;  see  Sec.  IV-A  for  further 
discussion.  For  beam  energies  of  several  MeV,  =  1  for  t^  expressed  in 
equivalent  mils  of  titanium. 

B.  Foil  Focusing 

A  foil  in  the  path  of  an  intense  electron  beam  quickly  becomes 
conducting,  thereby  shorting  the  electrostatic  repulsive  field  of  the 
beam.^®~^^  This  momentary  shorting  imparts  a  focusing  impulse  much  like  a 
thin  lens.  For  paraxial  beams,  the  foil  focal  length  varies  as^^ 
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(24) 


yIa 

“  ®f®rms  (l-fi)Ib  ’ 

where  a  is  the  rms  beam  radius,  g,  -  0.4  is  a  geometrical  coefficient 
rms  I 

determined  by  the  beam  radial  profile  and  wall  radius,  and  f^  is  the  gas 
charge-neutralization  fraction  averaged  over  the  two  sides  of  the  foil. 

Foil  focusing  has  little  effect  if  it  is  weak, 

>  Xp/4  ,  (25) 

or  if  the  IFR  cell  is  short, 

»  z.  (26) 

Here  is  the  IFR  betatron  wavelength  given  in  Eq.  (9).  Condition  (26)  is 
always  met  if 

<27) 

If  none  of  these  conditions  is  met,  the  beam  is  over-focused  and  some 
electrons  are  pushed  to  outside  a^.  This  increases  the  effective  injection 
radius,  and  it  adds  beam  wings  to  the  radial  profile.  The  over-focusing 
can  also  produce  loss  of  transported  current  as  the  beam  scrapes  the  pipe 
walls.  All  of  these  effects  are  strongest  near  the  peak  of 

Another  adverse  foil  effect  is  that  the  focal  length  f^  varies  with 
particle  position.  This  enharmonic  focusing  raises  the  beam  emittance  by^<^ 
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where  ~  0. 1-0.5  is  a  geometrical  coefficient  determined  by  the  beam 
profile  and  pipe  diameter.  Foil  focusing  thus  produces  inverse  emittance 
tailoring  during  the  rising  portion  of  the  beam  current.  Such  inverse 
tailoring  must  be  eliminated  to  produce  radius  tailoring.  The  simplest 
solution  is  to  make  the  entrance  foil  sufficiently  thick  that  emittance 
growth  from  scattering  dominates  the  growth  from  anharmonic  focusing.  This 
requires,  according  to  Eqs.  (23)  and  (28),  a  foil  thickness  satisfying 

kjtj  »  laj(l-f.)Ijj/I^I^  .  (29) 

One  can  show  that  this  is  compatible  with  the  emittance  required  in  Eq. 

(20)  provided 

.  (30) 

In  general,  foil  focusing  can  be  disregarded  if  the  beam  is  highly 
paraxial,  «  yI^* 

C.  Space-Charge  Depression 

Foil  focusing  stems  from  charge  induced  on  the  foil .  This  charge  not 
only  shorts  the  radial  electric  field,  but  it  also  produces  an  axial 

43 

electric  field.  We  show  here  that  the  axial  field,  not  the  radial  field, 

is  responsible  for  space-charge  depression  of  the  beam  energy. 

As  a  beam  leaves  a  conducting  foil,  it  decelerates  in  response  to  the 

(image)  charge  on  the  foil.  The  induced  axial  field  is  electrostatic  and 

given  by  E^^  =  -  34»^/3z,  where  is  the  potential  from  the  foil  charge. 

The  field  E  ,  lowers  the  beam  energy  by 
zt 
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(31) 


2  7  r  3i.  Hf 

ACvmc  )  =  e  Jdz  Jdr  3^  ^  . 


0  0 


Here  ijj(r)  is  the  fraction  of  beam  current  within  radius  r,  and  z  =  0  is 
the  foil  location.  The  foil  potential  satisfies  Laplace's  equation 
outside  the  foil  and  is  given  on  the  foil  by  the  neutralization  requirement 


!!i 

ar 


=  E 


ro 


(32) 


where  is  the  radial  electric  field  from  the  beam  and  IFR  plasma. 


The  foil  fields  affect  the  beam  over  dz  ~  3a  outside  the  foil.^^  If 

rms 

f^  »  3a^|^g,  the  beam  profile  can  be  taken  as  constant  in  z,  and  the  order 
of  integration  in  Eq.  (31)  can  be  interchanged  to  yield 


A(Ymc^)  =  e  Jdr  ^  Jdz 


aT 


=  -  e  Jdr  4'f(r,0,t) 

0 

b 

=  e  dr  i.E 
J  b  ro 

0 

(l-f)Ib  •  (33) 

0 

Here  we  have  integrated  by  parts,  applied  boundary  condition  (32)  plus  the 
condition  that  =  0  at  r  =  b  and  z  =  ®,  imposed  ij^  =  0  at  r  =  0,  and 
expressed  E^^  in  terms  of  an  average  neutralization  fraction,  l-f(r).  For 


=  -  dr 


J- 


2ei, 


rc 


!5 


a  flat-topped  beam  in  vacuum,  result  (33)  reduces  to  Ay  = 

2  2 

-  lln(b  /a^  )+0.5]  (I^/I^).  One  can  show  that  this  is  just  the  energy 
needed  to  regenerate  the  beam  self-fields  (magnetic  and  electrostatic) 
outside  the  foil. 

Although  the  energy  loss  Ay  is  eventually  recovered  as  the  beam  passes 

through  the  IFR  exit  foil,  the  value  of  y  within  the  IFR  cell  remains 

reduced.  This  alters  the  radius  tailoring  because  y  varies  with  time  t 

(and  with  radius  r  because  varies  with  r).  These  effects  become 

pronounced  for  Ay  -»  -(y^-l),  where  y^  is  the  injection  energy.  This  value 

43 

of  Ay  defines  the  space-charge  limiting  current  in  the  IFR  cell.  Above 
this  current,  the  beam  slows  and  expands  to  fill  the  pipe,  until  scrape-off 
or  low  axial  velocity  reduces 

IV.  PARAHETER  VARIATIO^^:  IN  A  PASSIVE  IFR  CELL 


A.  Design  Equations 

In  this  section,  we  use  the  previous  results  to  construct  a  set  of 

algebraic  equations  for  designing  a  passive  IFR  tailoring  cell  for  a  given 

beam.  For  simplicity,  we  ignore  space-charge  depression  and  foil  focusing 

(valid  for  highly  paraxial  beams),  and  we  assume  that  the  beam  approaches 

the  IFR  cell  with  constant  energy  y,  emittance  ,  and  radius  a^.  The 

adjustable  IFR  parameters  are  the  entrance  and  exit  foil  thicknesses,  the 

cell  dimensions  b  and  2,  and  the  gas  composition  and  pressure  P.  The 

output  is  a  prescribed  variation  in  the  beam  radius  given  as  follows:  a^^ 

is  to  start  collapsing  at  time  Tj,  fall  until  time  and  saturate  at 

a  /■f2  for  T  >  T  .  The  time  t.  marks  the  transition  from  free  expansion  to 
o  ~  e  1 
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pinched  equilibrium,  and  the  time  is  where  emittance  growth  prevents 
further  beeun  collapse.  The  actual  transitions  at  and  x^  are,  of  course, 
gradual  rather  than  abrupt. 

The  output  beam  temperature  is  specified  in  terms  of  a  maximum  value, 

T  ,  the  choice  of  which  is  discussed  later.  This  bound  limits  the  IFR 
ndx 

pinch  strength  to 


'i^b  ^ 


max 


2cT  /e  , 
max  ’ 


(34) 


according  to  Eq.  (4).  Here  is  an  equivalent  pinch  "current" 

characterizing  the  strength  of  the  ion  charge  seen  by  the  beam.  Because 

f.I.  increases  monotonically  with  time,  it  reaches  I  at  the  end  of  the 
ID  max 

beam  pulse,  x  *  Xp.  Applying  this  bound  to  Eq.  (21)  yields 
'’i  ■ 

This  equation  specifies  the  gas  ionization  rate  v^,  and  thus  the  gas 

pressure  P,  in  terms  of  the  maximum  ion  pinch  "current"  I  and  the  beam 

max 

charge  at  the  times  x^,  x^,  and  Xp. 

Ve  now  compute  the  thickness  of  the  IFR  entrance  foil.  The  beam  radius 
at  time  x^  equals  a^/{2.  To  produce  this  radius,  Eqs.  (15b)  and  (16) 
indicate  that  the  beam  emittance  should  satisfy 

Cq  =  jYVilQb(x^)  ^  Qj,(x.)]/8I^  .  (36) 

According  to  Eq.  (23),  the  entrsmce  foil  thickness  should  therefore  satisfy 

kjtj  =  YVj[Qb(Xg)+Qb(x.)|/4I^  -  2(E./a^)2  ,  (37a) 
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where  the  original  beam  emittance  is  assumed  to  be  small,  e.  <  Ve  can 

rewrite  condition  (37a)  using  Eq.  (35)  as 


kftf 


yI 

max 

r  Qb<%> "  Qb<^i>  1 

o 

fe. 

1 

K<^p>  -  “b<^e)  * 

—  z 

a 

L  oj 

(37b) 


The  length  z  o£  the  IFR  cell  determines  the  amount  of  radius  tailoring 
achieved  and  its  onset  time  The  length  is  given  from  Eqs.  (16)  and 

(35)  by 


r  Yl, 


z  = 


.1/2 


max 


fQ^^(Tp)  -  0.5  Qj^(t^)  +  0.5  Qj,(T.)'l 


V^i> 


1/2 


(38) 


The  use  of  larger  z  would  lower  x^  and  give  more  tailoring.  We  point  out 
that  result  (38)  is  consistent  with  equilibrium  requirement  (9)  at  the 
tailoring-end  time  x^  only  if 


a  /a  .  >  n  , 

max  min  -  ’ 


(39) 


where  the  minimum  output  radius  a  .  =  a  /42. 

min  o 

The  maximum  head-to-tail  radius  variation  is  given  by 


a 

max 

a  . 
min 


Qb(%)  +  Qb<^i> 


,<"i> 


(40) 


If  desired,  one  could  specify  a  /a„.  as  a  design  objective  in  place  of 
X. ,  and  use  Eq.  (40)  to  calculate  Qj,(y.).  In  any  event,  the  pipe  radius  b 
is  chosen  to  minimize  beam  scrape-off: 


>  a  =  a 
max  o 


Qb<%)  "  V''i> 


2Qj^(Xi) 


(41) 
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Pipe  scrape-off  could,  of  course,  be  used  to  control  the  head  radius 

and  tailoring-onset  time  t^,  but  only  at  the  loss  of  current  in  the  head. 

Equations  (35),  (37),  (38),  and  (40)  are  a  prescription  for  designing 

an  IFR  radius  tailoring  cell.  To  convert  the  radius  tailoring  to  emittance 

tailoring,  the  beam  must  be  further  heated  at  the  IFR  exit  to  a  constant 

temperature,  T,  >  T  .  Emittance  tailoring  is  maximized  when  T,  »  T 

f  ~  max  *  f  max 

so  that  Ejj(‘t)  varies  directly  with  a|j(T).  According  to  Eqs.  (23)  and  (34), 
this  can  be  achieved  by  passing  the  beam  through  a  scattering  foil  whose 
thickness  t^  satisfies 

kjtj  .  2yrj/.c2  »  .  (42) 

The  beam  could  alternatively  be  heated  using  an  active  IFR  cell  or  similar 
enharmonic  centering  scheme.  The  IFR  exit  foil  should  then  be  thin. 

The  final  emittance  desired  for  the  beam  determines  T^.  To  keep  the 
radius  of  the  beam  body  unchanged  in  the  propagation  region  following  the 
IFR  cell,  we  set 

Tf  =  eleff/2c  ,  (43) 

where  is  the  effective  body  pinch  current  in  the  propagation  region. 

For  a  beam  propagating  in  dense  gas,  ^eff  equals  the  beam  current  minus  the 
plasma  return  current  enclosed  by  the  beam;  in  this  case,  Ig££  typically^^ 
lies  between  half  to  one  times  the  peak  beam  current  I^.  Choice  (43) 
matches  the  beam  body  but  not  the  head.  As  a  consequence,  the  overheated 
head  will  expand  as  the  beam  propagates.  The  head  mismatch  can  be  reduced 
by  raising  the  IFR  gas  pressure  P,  but  this  decreases  the  amount  of 
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emittance  tailoring  achieved.  At  times  T  >  the  emittance  is  always 
slightly  inverse  tailored. 

Before  closing  this  section,  we  point  out  that  the  design  formalism  is 
based  on  the  assumption  that  the  beam  should  taper  to  its  smallest  possible 
radius,  a^/JT.  In  some  instances,  however,  a  larger  minimum  may  be 
desired,  e.g.,  to  propagate  a  beam  stably  over  long  distances  outside  the 
IFR  cell.  In  that  case,  design  equations  (35),  (38),  and  (41)  can  be 


retained  by  setting  equal  to  the  pulse  duration  Xp.  But  Eq.  (36)  for 

1/2 

the  entrance  emittance  should  be  replaced  with  c  =  a  .  (yI  /2I.)  , 

^  n  min'  max  A'  ’ 


where  a  , 
min 


is  the  desired  minimum  beam  radius  at  time  x  . 

P 


B.  Examples  and  Comparison  with  Simulation 

In  this  section  we  illustrate  the  radius  variation  produced  by  a 
passive  IFR  cell.  The  beam  is  assumed  to  have  the  half-sine  current  form 
given  by  Eq.  (12a),  with  nominal  parameters  I^  =  17  kA,  *  1  cm,  c.  =  0.3 
rad-cm,  *  10,  and  Xp  *  40  nsec.  Let  us  design  the  IFR  cell  to  produce  a 
tailoring  onset  time  x^  =  4  ns,  a  tailoring  end  time  x^  =  Xp  =  40  ns,  and  a 
peak  ion-channel  strength  I^^^  =  2.2  kA.  Using  the  fact  that  Qjj(t)  = 
(I^Xp/n) ( l-cos(nx/Xp) 1 ,  we  compute  the  necessary  gas  ionization  rate  from 
Eq.  (35)  as  *  10^  s~^,  corresponding  to  a  gas  pressure  P  =  10  mtorr  in 
air.  Equation  (37b)  indicates  that  the  IFR  entrance  foil  should  be  0.5  mil 
of  titanium,  corresponding  to  a  net  injection  emittance  of  =  0.6  rad-cm. 
Equation  (38)  indicates  the  cell  should  be  40  cm  long,  and  scrape  condition 
(41)  indicates  that  the  wall  radius  b  should  exceed  5  cm.  We  set  b  =  10  cm 
for  operational  flexibility.  These  beam  and  IFR  cell  parameters  are 

9 

similar  to  those  used  in  the  SuperlBEX  experiments. 
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The  predicted  beam  radius  at  cell  exit  is  given  by  Eq.  (15b), 

subject  to  the  collapse  restriction,  a^^  >  a^/-(2.  Although  foil  focusing  is 

predicted  to  be  strong  for  this  beam,  we  ignore  it  here.  The  effect  is 

incorporated,  however,  in  the  simulations  described  in  Ref.  37. 

Figure  1  plots  a^^CT)  from  Eq.  (15b)  for  several  values  of  the  IFR  gas 

pressure  F.  The  pressure  P  affects  both  the  tailoring  onset  time  and 

the  end  time  T  ,  but  does  not  affect  the  head  radius  a  .  The  latter  is 
e  max 

limited  by  free  expansion  to  a^^^^  =  4.8  cm  in  this  case.  Raising  P 

shortens  T.  and  so  that  the  beam  pinches  more  quickly.  In  the  beam  tail 

(t  >  30  ns),  the  radius  is  roughly  constamt  because  the  beam  current 

decreases  rapidly  and  little  additional  ionization  is  produced. 

Figure  2  shows  the  effect  of  changing  the  input  emittance  (the 

entrance  foil  thickness)  while  holding  P  =  10  mtorr.  The  emittance 

affects  the  tailoring  end  time  x  and  the  head  radius  a„^„,  but  does  not 

e  max 

affect  the  tailoring  onset  time  Xj.  Little  tailoring  occurs  if  is  too 

low,  and  the  beam  does  not  pinch  tightly  if  is  too  high. 

Figure  3  shows  the  effect  of  the  cell  length  z  on  tailoring.  The  cell 

length  affects  both  a  and  x.  but  usually  has  little  effect  on  x^. 

^  max  1  e 

Consequently,  raising  z  increases  tailoring  mainly  in  the  beam  head.  In 
general,  the  cell  should  be  made  too  long  rather  than  too  short  to  maximize 
tailoring  flexibility. 

In  Fig.  4  we  correct  for  space-charge  depression  of  the  beam  energy 
using  Eq.  (33).  The  reduction  in  y  approaches  50%  near  peak  current  in 
Fig.  4a  for  the  standard  case  described  above.  Nevertheless,  the  effect  on 
radius  tailoring  is  relatively  modest,  as  shown  by  the  comparison  in  Fig. 

4b  between  the  corrected  beam  radius  (solid  line)  and  the  uncorrected 
radius  (dashed  line).  However,  if  the  current  were  to  approach  the  space- 
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charge  limiting  value  ~  30  kA,  the  beam  would  expand  to  fill  the  pipe  with 
the  radius  virtually  untailored. 

To  test  our  model,  ve  have  compared  its  predictions  with  those  from  a 

37 

fully  electromagnetic  particle  code,  FRIEZR.  In  particle  codes  (and 

experiments),  the  beam  profile  rounds  so  that  the  half-current  radius  ^2/2 

or  rms  radius  a  is  a  more  useful  measure  than  the  edge  radius  a.  .  To 
rms  D 

allow  for  this  difference,  we  divide  the  analytical  result  for  aj^  by  'fl, 

which  is  the  appropriate  scale  factor  for  both  a. and  a  for  a  flat- 

topped  beam.  In  Fig.  5,  we  plot  the  analytical  radius  aj^/'fz  (dashed  line) 

and  the  simulation  radii,  a,.,  (solid  line)  and  a  (dot-dashed  line),  for 

i/z  rms 

a  beam  injected  with  a  Gaussian  radial  profile  and  the  nominal  parameters 
given  above.  The  analytical  result  tracks  ^^^2  closely,  especially 

for  T  >  10  ns.  The  discrepancy  at  early  times  occurs  because  the  beam  is 
not  flat-topped  and  because  it  over-expands  without  equilibrating. 

Overall,  the  agreement  between  code  and  theory  is  good,  indicating  the 
utility  as  well  as  the  limitations  of  the  model. 


V.  BEAM  CENTERING 


A.  Passive  IFR  Cell 

Another  major  benefit  of  IFR  propagation  is  suppression  of  azimuthal 
asymmetries  like  transverse  offsets  (m=l),  non-circular  profiles  (m=2),  and 
filamentation  (m  >  2).  In  all  cases,  the  price  of  symmetrization  is 
increased  beam  emittance.  Here  we  concentrate  on  beam  centering. 

There  are  two  types  of  centering  forces  from  an  IFR  cell:  those 
produced  by  the  IFR  channel  itself,  and  those  from  the  pipe  walls.  The 
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vail  forces  are  easily  computed  if  we  assume  that  the  pipe  is  perfectly 

conducting  and  that  the  ion  channel  is  cylindrical  and  long,  b  9/3z  <<  1. 

The  latter  condition  is  usually  well  satisfied  if  the  beam  is  paraxial  and 

the  IFR  cell  is  long,  z  »  b.  The  ion  channel  can  then  be  treated  as  a 

long  cylinder  of  line  charge  q..  If  the  channel  resides  at  a  location 

displaced  from  the  pipe  center,  it  induces  a  charge  on  the  pipe  that  can 

2 

be  represented  as  an  image  line  charge  of  -q^  located  at  (b/y^)  5^^  outside 
52 

the  pipe.  This  image  charge  deflects  the  beam  with  an  average  force 
given  by 


~y  " 


2eq. 


((b/yc)V  i^b)’ 


(44) 


where  is  the  beam  centroid.  Deflection  forces  from  the  beam  image 

_2 

charge  and  current  cancel  to  order  y  • 

In  a  passive  IFR  cell,  the  channel  is  essentially  centered  on  the  beam, 
with  a  charge  q.  =  v^Qj^/c.  Equation  (44)  can  then  be  written  as 


F  =  - 

-y 


2eN>iQb 


(b^ 


2, 

yK  )c 


^b  • 


(A5) 


2 

For  yj^  <<  b,  the  centering  force  reduces  to  F^  =  (2ev^0j^/b  c)y|^.  This 
indicates  that  the  beam  centers  in  a  characteristic  distance 


z^  =  (n/2)  'IyI^/2v.Qjj  b  .  (46) 

Phase  mixing  in  the  anharmonic  pinch  force  from  q^  damps  subsequent 
oscillations  in 

Result  (46)  neglects  any  separation  between  the  beam  and  ion  channel. 
Such  separation  will  in  fact  develop  because  the  rise  in  with  t  causes 


the  beam  head  to  center  more  slowly  than  the  body.  This  difference  tilts 
the  beam  and  causes  the  plasma  channel  to  lag  behind  and  separate,  thereby 
generating  additional  beam  deflection  forces.  Nevertheless,  Eq.  (46) 
provides  a  useful  estimate  for  the  minimum  cell  length  needed  to  center  a 
beam  in  a  passive  IFR  cell. 

IFR  centering  has  been  demonstrated  using  the  SuperlBEX  beam  described 

earlier.  The  centering  experiments  were  conducted  in  a  pipe  114  cm  long  by 
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4.7  cm  in  diameter  and  filled  with  20  mtorr  of  nitrogen.  Nearly  all  of 
the  beam  was  centered  by  the  pipe,  as  predicted  by  Eq.  (46).  The  weakly 
pinched  head,  however,  was  lost  to  pipe  scrape-off  as  expected. 

B.  Centering  in  an  Active  IFR  Cell 

In  an  active  IFR  cell,  the  channel  is  created  externally  and  usually  in 
the  center  of  the  pipe.  A  centered  channel  eliminates  the  wall  force,  as 
can  be  shown  by  setting  =  0  in  Eq.  (44).  This  leaves  the  channel  to 
guide  as  well  as  pinch  the  beam. 

Unlike  wall  centering,  channel  centering  cannot  be  computed  simply  by 
evaluating  the  static,  net  deflection  force  on  the  beam.  Instead, 
dynamical  forces  from  orbit  phase-mixing  and  beam  distortion  must  be 
included  as  well.  To  illustrate  this,  consider  a  hollow  beam  surrounding  a 
narrow  channel.  Because  the  beam  produces  no  fields  or  forces  within 
itself,  there  is  no  force  on  the  enclosed  channel,  and  thus  there  is  no  net 
deflection  force  on  the  beam  by  reciprocity.  However,  the  beam  does  not 
maintain  its  original  structure  but  collapses  towards  the  channel.  This 
collapse  causes  the  beam  to  fill  in  and  distort  until  it  forms  a  new 
profile  centered  about  the  channel.  This  is  a  dynamical  process  that 
cannot  be  predicted  using  a  static,  rigid-rod  force  calculation. 
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In  fact,  narrow  channels  are  especially  effective  at  centering  beams, 

hollow  or  not.  The  main  reason  is  that  the  channel  force  is  then  highly 

enharmonic,  causing  the  beam  electrons  to  orbit  about  the  channel  but  not 

in  phase  with  one  another.  As  a  consequence,  the  orbits  quickly  lose 

coherence  and  phase-mix  such  that  the  non-rigid  beam  distorts  into  a  new 

profile  centered  about  the  channel.  Phase-mix  centering  can  occur  in 

distances  as  short  as  a  few  X-  for  a  highly  anharmonic  force.  A  second 

P 

reason  narrow  channels  are  effective  at  beam  centering  is  that  the  full  ion 
charge  attracts  each  beam  electron,  thereby  maximizing  the  channel  force. 

To  illustrate  beam  centering  in  an  active  IFR  cell,  we  consider  a  broad 
beam  and  a  wire-like  ion  channel  created  externally  with  a  fixed  line 
charge  q^  per  unit  length.  The  line  charge  q.  attracts  each  beam  electron 
with  a  force 

Fj.(r)  =  -  2eq./r  ,  (47) 

where  r  is  the  radial  distance  from  the  channel.  This  force  heats  the 
electrons  until  pressure  balance  is  achieved,  3(nj^Tj^)/3r  =  nj^F^.,  where  n^^ 
is  the  beam  density  and  Tj^  is  the  transverse  beam  temperature.  Using 
pressure  balance  and  Eq.  (47),  the  average  final  temperature  is  given  by 

b  b 

■  "b"'  -  -  N,,''  Jdr  (n^Tj^) 

0  0 

b  b 

=  -  Jdr  nr^  n^^F^  =  eq.N^^"^  Jdr  2nrn^^ 

0  0 

=  eq^  ,  (48) 
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All  electrons  thus  gyrate  about  the  channel,  but  they  gyrate  at  different 

frequencies  determined  by  their  initial  location  r  from  the  channel. 

The  variation  in  oscillation  time  causes  the  electrons  to  become  out  of 

phase  with  one  another.  As  the  particles  lose  coherence,  the  beam  evolves 

to  a  smooth,  stable  distribution  centered  about  the  channel.  That  is, 

phase  mixing  of  the  particle  orbits  damps  fluctuations  in  the  beam  profile 

12 

and  centroid  location.  Without  this  damping,  the  fluctuations  would 
persist,  and  the  beam  would  forever  oscillate  about  the  channel.  The 
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latter  scenario  occurs  if  the  channel  force  is  perfectly  harmonic,  F^(r)  “ 
r,  so  that  t^  is  independent  of  r.  For  this  reason,  a  flat-topped  ion 
channel  broader  than  the  beam  is  ineffective  at  beam  centering. 

Electrons  starting  far  from  the  channel  take  longer  to  cross  it,  and 
thus  they  phase-mix  more  slowly  than  electrons  starting  near  the  channel. 
The  characteristic  distance  for  phase-mix  centering  is  therefore  governed 
by  the  longest  betatron  wavelength,  given  from  Eq.  (50b)  by 

'  2(a^.y^)Jl.rIj/cqj  ,  (51) 

where  a^  is  the  beam  edge  radius  and  yj^  is  the  beam  offset  at  injection. 

The  beam  first  crosses  the  channel  at  z  =  but  several  are 

c  P  ’  3 

needed  for  phase-mixing  to  damp  subsequent  oscillations  of  the  centroid. 

If  the  electrons  are  injected  warm  rather  than  cold,  their  tangential 
velocity  causes  the  orbits  to  process  about  the  channel.  This  precession 
produces  a  beam  that  is  azimuthally  symmetric  as  well  as  centered. 
Asymmetries  like  beam  filamentation  and  non-circularity  are  thereby 
removed . 

Symmetrization  and  centering  come  at  the  expense  of  beam  emittance.  A 
lower  bound  for  the  emittance  can  be  derived  by  combining  temperature 
result  (48)  with  an  estimate  for  the  final  beam  area,  yielding 


min 

e 

n 


(52) 


This  is  consistent  with  Eqs.  (17)  and  (18)  for  yjj  =  0  and  f^Ij^  =  2cqj. 
Note  that  the  offset  is  unimportant  unless  it  is  large,  y^^  >  a^. 

An  active  IFR  cell  could  be  combined  with  a  separate  passive  IFR  cell 
to  produce  a  beam  that  is  both  centered  and  emittance  tailored.  The  most 
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efficient  arrangetaent  is  to  place  the  active  IFR  cell  last,  so  that  the 
heating  from  this  cell  eliminates  the  need  for  a  final,  thick  scattering 
foil.  In  the  reverse  order,  the  emittance  gained  from  the  centering  cell 
cannot  be  efficiently  utilized. 

A  limitation  of  active  IFR  cells  is  that  they  do  not  become  fully 
operational  until  the  beam  current  rises  above  cq^.  In  the  early  beam 
head,  this  condition  is  not  met,  and  the  head  propagates  in  an  overdense 
plasma.  This  region  is  often  sufficiently  short,  however,  that  there  is 
little  adverse  impact  on  the  beam. 


VI.  PLASMA  MOTION 


A.  Ion  Channel  Contraction  and  Expansion 

The  assumption  that  the  ions  are  stationary  and  the  plasma  electrons 
are  instantaneously  ejected  is  not,  of  course,  strictly  justified.  In  this 
and  the  following  sections,  we  discuss  ion  pinching,  ion  hose  motion,  ion 
impact  ionization,  and  magnetic  trapping  of  the  plasma  electrons. 

The  beam  electrostatically  pinches  the  ion  channel,  causing  it  to 
contract.  Such  contraction  takes  place  on  an  ion  crossing  time  given  by 

=  {rta^/2c)(M/m)"^^(I^/Ijj)^^^  ,  (53) 

where  M  is  the  ion  mass.  Contraction  is  unimportant  if  is  much  larger 
than  the  beam  pulse  duration  t^,  or  if  the  channel  is  narrow  as  was  assumed 
in  the  analysis  of  an  active  IFR  cell.  If  the  channel  is  broad,  its 
contraction  raises  the  ion  charge  f^Ij^/c  seen  by  the  beam;  this  raises  the 
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equilibrium  value  of  as  indicated  by  Eqs.  (3)  and  (A).  Using  the 

arguments  given  for  the  beam,  channel  contraction  is  ultimately  limited  to 
-1/42,  and  a  broad  channel  that  contracts  by  this  amount  doubles  fjljj  and 
causes  to  increase  by  tT?.  However,  the  contraction  can  be  larger  at 

X  =  x^,  depending  on  how  enharmonic  the  channel  pinch  force  is. 

The  mai».  effect  of  channel  contraction  in  a  passive  IFR  cell  is  to 

raise  e  if  the  tailoring-end  time  t  <  t  ,  or  to  reduce  x  if  t  >  t  .  In 
n  e  c  e  e  c 

either  case,  the  effect  is  usually  modest  because  the  channel  is  being 
continuously  created  by  the  beam.  Channel  collapse  is  therefore  generally 
unimportant  in  passive  IFR  cells,  as  well  as  in  narrow-channel  active  IFR 
cells. 

Late  in  the  beam  pulse,  the  beam  current  drops  and  the  beam  density 
ultimately  falls  below  the  ion  density,  f.  >  1.  Thermal  pressure  and 
mutual  repulsion  then  drive  the  ion  channel  to  expand.  Such  expansion 
could  cause  the  tail  of  the  beam  to  flare,  especially  if  the  beam  current 
decays  over  a  time  longer  than  an  ion  crossing  time  In  a  typical 

passive  IFR  cell  where  fj  «  1  in  the  beam  body,  flaring  is  not  expected 
until  very  late  in  the  pulse.  Tail  flaring  is  similarly  unimportant  in 
active  IFR  cells  provided  the  ion  channel  is  initially  much  smaller  than 
the  beam.  Expansion  and  contraction  of  the  ion  channel  can,  of  course,  be 
reduced  by  using  a  fill  gas  with  a  large  ion  mass  M. 


B.  Ion  Hose  Instability 

Potentially  more  disruptive  than  channel  contraction  and  expansion  is 

j  u  16,20,21,38-40  . 

transverse  motion  caused  by  the  lon-hose  instability.  This 

is  an  electrostatic  instability  produced  as  the  beam  and  ion  channel  move 
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sideways  and  separate.  The  instability  convects  back  from  the  beam  head 


in  a  characteristic  growth  time  At  ~  over  a  growth  distance  tz  ~  \^/2n. 

The  instability  tends  to  saturate  as  the  beam  offset  becomes  comparable  to 

the  beam  and  channel  radii.  The  offset  is  then  converted  to  beam  emittance 

through  phase-mixing  of  the  particle  orbits.  The  net  effect  is  a  centered 
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beam  but  with  a  flared  radius  in  the  tail  (assuming  there  is  no 
downstream  re-excitation  of  the  instability  and  >  2Tp/T^).  Tail 

O 

flaring  was  in  fact  observed  on  ATA. 

Because  an  active  IFR  cell  is  typically  several  long,  tail  flaring 
from  ion  hose  can  be  appreciable  unless  the  pulse  length  is  short,  x^  <  x^ 
for  initial  offsets  *-0.1  a^^.  Ion  hose  can  be  suppressed,  however,  by 
making  the  channel  strength  rise  with  time  x;  for  example,  the 

ionizing  laser  pulse  can  be  made  to  overlap  the  beam  pulse.  A  rise  in 
detunes  the  instability  and,  according  to  recent  analysis  to  be  reported  in 
a  later  paper,  permits  pulse  lengths  up  to  5  x^  or  longer. 

In  a  passive  IFR  cell,  the  cell  length  is  typically  <  X^,  but  this  is 
still  sufficient  for  substantial  hose  growth.  In  this  case,  however,  the 
channel  s*  ^ngth  inherently  rises  with  x,  and  the  channel  also  follows  the 
beam  throug.  beam  ionization.  Pulse  lengths  of  5  x^  or  longer  are  thus 
permitted  before  ion  hose  and  tail  flaring  become  excessive. 

C.  Ion-Impact  Ionization 

If  the  ions  become  sufficiently  energetic,  they  will  ionize  the  gas. 
This  additional  ionization  increases  f.,  thereby  altering  the  beam  pinch 
force  and  possibly  threatening  the  requirement  that  f^  <  1  in  the  beam 
body.  Hence,  ion-impact  ionization  should  in  general  be  minimized  for  IFR 
cells  to  function  properly. 

Ion  ionization  competes  with  beam  ionization  if  n^v^^  >  where 

is  the  ion  ionization  rate.  This  can  be  rewritten  as 
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(54) 


where  a.  is  the  beam  ionization  cross  section,  a,  is  the  ion  ionization 
D  1 

cross  section,  and  is  an  average  ion  speed. 

The  beam  ionization  cross  section  is  given  from  the  Born-Bethe  model  at 
relativistic  energies  by^^ 


«Tb(Y)  =  4iir^  [1  +  Cj  ln(Y)]. 


where  r  ■  3.86x10”  cm  is  the  electron  Compton  wavelength,  is  the 
o  o 

nuclear  charge  of  the  target  molecule,  and  >  1  and  e  0.2  are  gas 

constants.  The  ion  ionization  cross  section  can  be  approximated  using 
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Pirsov's  model  by 


aj(Vj)  =  6.3xl0”^^(Z.+Z^)(v./c)(Rjj/Wj)  cm^ 


(56a) 


in  the  velocity  range  5  <  (v./v^)  <  50,  Here  Z^  is  the  nuclear  charge  of 
the  incident  ion,  is  the  ionization  energy,  Rj^  =  13.6  eV  is  the  Rydberg 
energy ,  and 


V  =  5.4x10  ^(Z.-.-Z„)”^^^(W./R.  )c 

o  10  in 


(56b) 


is  a  scale  velocity.  Substituting  approximations  (55)  and  (56a)  into 
condition  (54)  indicates  that  ion  ionization  becomes  important  when 


c  z 

f.  >  3x10"^  (V./Rjj)(c/v.)^ll  +  Cj  ln(Y)]  • 

1  -  ^i"*"  o 
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Energy  conservation  restricts  the  ion  energy  to  the  available 


electrostatic  field  energy: 


< 


e(l-fi)Ib 


(58) 


vhere  r  <  a^^  is  the  starting  position  of  the  ion.  Charge  exchange  or 
similar  processes  can  reduce^^  the  energy  to  below  limit  (58),  but  this 
cooling  is  usually  unimportant  at  the  gas  densities  used  in  IFR  cells. 
Averaging  the  ion  energy  over  both  time  and  over  all  starting  positions 
yields  a  mean-squared  speed  given  by 

<VjS  <  e(l-fj)Ijj/2Mc  ,  (59) 


for  an  ion  channel  initially  uniform  out  to  radius  aj^. 

Inserting  limit  (59)  into  condition  (57)  indicates  that  ion  ionization 
is  unimportant  relative  to  beam  ionization  if 


c  z 

(i-fj)fjijj  <  6x10"^  rTz~  ^1 

i  o 


=  0.5  Z  (W,/R.  )[1  +  0.2  ln(Y)I  kA  . 
o  in 


(60a) 


Here  we  have  set  the  ion  mass  M  equal  to  (Z^^+Z^)  times  the  proton  mass, 
assumed  Zj  =  Z^,  and  have  set  =  2.5  and  =  0.2.  Because  Firsov's 
model  often  underestimates^^  <t.  by  a  factor  of  two,  we  shall  impose  the 
more  restrictive  requirement 


(l-f.)f.Ijj  <  (Z^/4)(W./Rj^)(l+0.2  ln(Y)]  kA  .  (60b) 

This  condition  is  always  met  if  either  <  Z^  kA  or  <  Z^/4  kA,  for  an 

ionization  energy  W.  =  13.6  eV. 
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Condition  (60b)  indicates  that  the  neglect  of  ion  ionization  is 
justified  provided  is  large  and  the  beam  current  or  channel  strength  are 

modest.  For  example,  in  argon  (Z^  °  18,  =  15.8  eV),  ion  ionization  can 

be  ignored  if  ^  20  kA  or  ^  ItA;  using  a  heavier  gas  would  raise 

these  limits,  as  well  as  reduce  other  adverse  effects  from  ion  motion. 

Even  when  ion  ionization  is  important,  the  increase  in  f^^  is  usually  not 
the  explosive  growth  described  by  Olson, because  the  secondary  ions  are 
created  near  the  beam  axis  where  the  primary  ion  energy  is  high  but  the 
electric  field  is  low.  Consequently,  the  secondary  ions  gain  little  field 
energy  and  thus  contribute  little  further  ionization.  Note  that  ion  and 
beam  ionization  can  be  eliminated  altogether  in  an  active  IFR  cell  by 
operating  at  gas  densities  sufficiently  low  that  )'fp 

Condition  (60b)  is  then  not  needed. 

D.  Trapping  of  Plasma  Electrons 

In  modeling  the  IFR,  we  have  assumed  that  the  plasma  electrons  are 
instantaneously  expelled  to  outside  the  beam.  This  is  usually  well 
justified  for  intense  beams,  unless  the  beam  current  is  so  high  that  it 
magnetically  traps  the  plasma  electrons.  In  this  section,  we  examine 
magnetic  trapping  and  the  bounds  it  imposes  on  IFR  propagation.  A  similar 

4 

analysis  was  given  by  Briggs  and  Yu. 

Consider  for  simplicity  a  long  beam  and  ion  channel  with  identical 
radial  profiles.  The  beam  and  channel  then  produce  a  radial  electric  field 
given  by 

z 

and  an  azimuthal  magnetic  field  given  by 


(62) 


Be(r)  -  ^  K(r)  . 


rc  "b' 


Here  is  the  axial  beam  velocity,  and  ijj(r)  is  the  fraction  of  beam 
current  within  radius  r. 

The  fields  and  Bg  alter  the  energy  and  axial  momentum  of  a  plasma 
electron  according  to 


d  /  2.  ^  „ 

dt  ^ 


(63) 


and 


ff  <  ■  ®St®e  • 


(64) 


where  y  is  the  mass  factor,  6c  is  the  axial  velocity,  and  6  c  is  the 

F  pz  pr 

radial  velocity  of  the  electron.  The  total  time  derivative  satisfies 


ft  o  i- 

dt  “  Ppr®  3r  * 


(65) 


Solving  these  equations,  assuming  constant  f^  and  and  zero  energy 
initially,  yields 


Yp  =  1  .  (l-f,)LIj,/6,I^ 


(66) 


and 


y„6„,  =  LI. /I.  , 

p  pz  b  A  ’ 


(67) 


where  I.  =  mc'^/e  and 
A 


r' 


L(r,r')  «  2  i^(r") 


(68) 
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is  a  dimensionless  parameter  related  to  the  distributed  inductance  and 
capacitance  of  the  beam.  Here  r  is  the  initial  position  of  the  plasma 
electron  and  r'  the  final  position. 

The  Lorentz  condition  restricts  the  plasma  mass  factor  to 


(1-0  ^-0  >  (1-0  2)-l/2 

'  pr  ^pz^  -  '  ^pz'' 


(69) 


Applying  this  bound  to  Eqs.  (66)  and  (67)  produces 


(1-f .)0 
'  z 


l-0.5f.-0.5/f.Y 
1  1  z 


<  1 


(70) 


2  2  -1/2 
where  y  f,  »  1  is  assumed  implicitly  and  y  =  (l-0_  )  is  the  beam 
Z  1  z  z 

axial  Lorentz  factor.  A  sufficient  condition  to  expel  all  plasma  electrons 
to  the  pipe  wall  at  r'  «  b  is  thus 


^i^b  -  2  ln(b/ajj)+l  ’ 


(71) 


where  L  has  been  set  to  its  maximum  value  for  a  flat-topped  beam.  In 
practice,  electrons  expelled  to  several  times  a^^  spend  most  of  their  time 
outside  the  beam,  and  they  therefore  have  little  effect  on  the  beam. 
Consequently,  magnetic  trapping  can  usually  be  ignored  for  ion  channel 
strengths  fjljj  <  5  kA. 


VII.  CONCLUSIONS 

In  this  paper,  we  have  modeled  IFR  cells  and  examined  their  utility  for 
conditioning  intense,  relativistic  electron  beams.  The  purpose  of  such 
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conditioning  is  to  prepare  the  beam  for  propagation  through  dense  gas 
without  disruption  from  the  resistive  hose  instability. 

Two  major  conditioning  goals  were  identified:  centering  and  emittance 
tailoring.  Centering  suppresses  the  hose  instability  by  reducing  the 
initial  beam  displacements  that  seed  the  instability,  while  tailoring 
reduces  the  hose  growth  rate  by  detuning  the  instability.  Ue  have  found 
that  wire-like,  active  IFR  cells  are  effective  at  centering,  and  that 
passive  IFR  cells  are  effective  at  tailoring.  Although  both  cells  have 
been  employed  for  these  purposes,  the  results  were  often  unsatisfactory  and 
not  well  understood.  Here  we  have  attempted  to  provide  a  simple  yet 
complete  design  procedure,  and  to  identify  the  major  limitations  of  the 
technique.  The  limitations  include  beam  focusing  from  the  entrance  and 
exit  foils,  magnetic  trapping  of  the  plasma  electrons,  and  ion  motion  in 
the  form  of  channel  collapse,  ion  hose  instability,  and  ion  ionization. 

We  conclude  from  our  analysis  that  IFR  cells  are  simple  yet  robust 

devices  that,  in  a  suitable  parameter  regime,  can  produce  beams  that  are 

both  well  tailored  and  centered.  As  has  already  been  demonstrated 
8  27 

experimentally,  ’  such  beams  propagate  much  mor«  stably  through  dense  gas 
than  unconditioned  beams. 
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Beam  radius  a.  as  a  function  of  x  for  three  differen 
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Fig.  2,  Beam  radius  a.  for  three  different  entrance  emittances:  t  =  0.3, 

D  o  ’ 

0.6,  1.2  rad-cm. 
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Fig.  4a.  Beam  current  and  space-charge  corrected  energy  y  as  functions 

of  T. 


Fig.  4b.  Beam  radius  ajj  with  space-charge  correction  (solid  line)  and 
without  space-charge  correction  (dashed  line). 
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Fig.  5.  Half-current  radius  3^/2  (solid  line)  and  rms  radius  (dot- 

dashed  line)  from  FRIEZR,  as  contrasted  with  the  analytical  result  a. /-fl. 
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APPENDIX;  DETUNING  THE  RESISTIVE  HOSE  INSTABILITY 


The  resistive  hose  instability^^  is  the  primary  factor  limiting  the 

propagation  of  intense,  self-pinched  electron  beams  in  dense,  un-ionized 

gases.  Resistive  hose  is  a  macroscopic  magnetic  kink  mode  that,  like  the 

ion-hose  instability  in  low-density  gas,  grows  and  convects  back  from  the 

beam  head  at  a  rate  determined  by  the  response  properties  of  the  beam  and 

plasma.  In  dense  gas,  the  plasma  response  is  characterized  by  a  magnetic 

diffusion  time  proportional  to  the  plasma  conductance  enclosed  by  the 

beam.  The  beam  response  is  characterized  by  the  betatron  wavenumber,  k„  = 

p 

2ll/Xp,  of  the  beam  electrons.  In  equilibrium. 


(Al) 


where  c  is  the  normalized  beam  emittance  and  I  re  <  II  is  the  effective 
n  erf  b 

beam  pinch  current.  Observe  that  there  is  no  explicit  dependence  on  beam 

radius  aj^  or  energy  t  in  Eq.  (Al). 

A  convective  instability  grows  with  time  t  into  the  pulse,  but  it 

eventually  decays  at  fixed  T  (provided  there  is  no  re-excitation  downstream 

of  the  accelerator).  As  a  consequence,  the  hose  amplitude  is  bounded  in  a 

pulse  of  finite  length.  Furthermore,  because  the  amplitude  scales  with  the 

initial  beam  displacement,  the  amplitude  can  be  reduced  by  placing  a 

8  26  29  36 

centering  cell  ahead  of  the  propagation  region.  ’  ’  ’  Centering  can  be 

achieved,  for  example,  using  the  active  IFR  cell  discussed  in  Sec.  V-B. 

A  second  means  of  reducing  hose  growth  is  to  "detune"  the  instability. 

8  30--33 

Detuning  is  possible  because  ’  the  growth  rate  of  the  instability 

depends  on  the  value  of  its  oscillation  wavenumber  k  relative  to  the  local 

value  of  ka.  Variations  in  k-  with  x  prevent  any  single  frequency  from 
P  P 
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remaining  in  resonance  throughout  the  pulse. 

Betatron  detuning  occurs  naturally  as  beam  ionization  causes  Ig££  to 

rise^^  from  zero  with  time  t.  According  to  Eq.  (Al),  this  detuning  can  be 

enhanced  by  tailoring  the  beam  emittance  to  fall  with  t.  The  amount  of 

tailoring  needed  is  dictated  by  the  number  of  hose  e-folds  allowed.  Recent 
28 

hose  theory  indicates  that  the  number  of  e-folds  varies  inversely  with  a 
detuning  coefficient, 


n 


eff  ^  a  ,  . 

I.  ^d  ax 


(A2) 


This  coefficient  is  usually  small  because  k.^  for  an  untailored  beam  rises 
on  a  monopole  decay  time  >>  x^.  An  untailored  beam  can  therefore  suffer 
many  orders  of  magnitude  of  hose  growth.  Raising  the  detuning  parameter  h 
by  two  or  more  through  emittance  tailoring  can  reduce  hose  growth 
dramatically. 

As  discussed  in  Sec.  IV,  a  beam  can  be  emittance  tailored  by  passing  it 

through  a  passive  IFR  cell  and  then  through  a  thick,  exit  foil  or  other 

heating  mechanism.  The  tailoring  profile,  and  thus  the  increase  in  q,  is 

limited  in  a  passive  IFR  cell  by  how  rapidly  the  beam  charge  rises  with 

X.  This  restricts  the  amount  of  hose  suppression  possible.  Other 

58  59 

techniques  proposed  for  tailoring  include  resistive  wire  cells,  ’ 
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current-carrying  wire  cells,  ’  ’  differential  focusing  cells,  and 

collapsing  ion  cells. 
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